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Low- lying excited states of the spin-1 antiferromagnetic Heisenberg 
chain with an impurity bond are investigated in terms of domain-wah ex- 
citations by means of analytical methods as well as a method of numerical 
diagonalization. It is shown that 1) the impurity bond brings about a mas- 
sive triplet mode in the Haldane gap, 2) the triplet state comprises three 
of the four ground states of an open chain, 3) the excited states with the 
domain wall trapped at the impurity bond, including the massive triplet 
mode, exhibit a Neel-type configuration of the magnetic moment around 
the impurity bond. On the basis of the present results, we propose a new 
energy versus magnetic-field diagram. 
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The spin-1 antiferromagnetic Heisenberg chain has been the subject 
of a large number of theoretical and experimental studies since Haldane's 
prediction^) of the difference between integer-spin and half-integer-spin an- 
tiferromagnetic Heisenberg chains. Most of the theoretical studies done so 
far have been concerned with the properties of the periodic chain. Recently 
Kennedy^) found that the open chain has a fourfold degenerate ground state 
composed of a singlet and a triplet which we call the Kennedy triplet, in 
contrast to a unique singlet ground state of the periodic chain. The four- 
fold degeneracy of the ground state, which was originally found in the 
so-called AKLT model^^ with open boundary conditions, is considered to 
reflect the hidden x symmetry in the open chain. 

In this letter we investigate low-lying excited states which include the 
Kennedy triplet of the spin-1 antiferromagnetic Heisenberg chain with an 
impurity bond. The impurity bond plays a mediating role between the 
open chain and the periodic chain. Using analytical methods as well as 
a method of numerical diagonalization, we calculate the energies of the 
low-lying excited states and also the magnetic moment at each site of the 
Kennedy triplet. In particular, exploring the dependence on the impurity- 
bond strength of the energy difference between the Kennedy triplet and 
the singlet is interesting and important, since it yields information on the 
transition of the ground state from the fourfold to the unique one. 

Before going into the details of the calculation, we briefly discuss our 
physical picture of the impurity-bond effect. It is convenient in many cases 
to express elementary excitations in the antiferromagnetic Heisenberg chain 
in terms of domain-wall excitations. If the excited state of the spin-1 pe- 
riodic chain is the sliding mode of domain-wall excitations, it is expected, 
analogously to the spin-1/2 Ising-like periodic chain, that the impurity 
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bond traps the domain wall and brings about the massive triplet mode be- 
low the bottom of the energy continuum. For the open chain, this massive 
mode should comprise three of the four ground states. It is also expected 
that the excited states with the trapped domain wall, including the mas- 
sive triplet mode, exhibit a Neel-type configuration of the magnetic moment 
around the impurity bond. 

We express the Hamiltonian which describes the spin-1 antiferromag- 
netic Heisenberg chain with an impurity bond as 

iV-l N 



hi,t = '^{S^S^, + S^Sy + A5|5|,) + 2K{Si ■ Si>f (A > 0) , (lb) 



where n represents the strength of the impurity bond, with At = and n=\ 
corresponding, respectively, to the open chain and the periodic chain. For 
convenience, we denote as the energy of the state corresponding to the 
singlet ground state at A = k = 1 and D = K = 0, and also as e['^^ the 
energies of the states corresponding to the Kennedy triplet at A = 1 and 
K= D = K = 0, where r = for the S^^^^^^ = Xl^Li*^! = state and r = ± 
for the -S^^otai = =1=1 states. Furthermore, we set Ao{k,) = ~ Es and 



We first calculate, by means of a variational method, Es and A^(k) 
(i = 0, 1) for the case of We introduce the following matrix-product- 

form wave functions $, \ and $i \ generalizing that of Kliimper 




(la) 



Ai{i^) = Ei^^-Es. 



et al.^^ which was originally proposed by Fannes et alJ^ The function $ 
describes the state with no domain wall and is defined as 
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where a^, C^, and (5^ are the spin states at the ^-th site, which correspond, 
respectively, to 51= 1, 0, and —1, and a± \={ax^i<Jy) / <Jxi Cy, and cr^ 
are the PauH matrices. On the other hand, the function $n describes the 
states with a domain wall at the n-th bond between the nth and (n+l)th 
sites, and is defined in terms of 0^ and an operator which may be called 
the wall operator, as 

= Trace [0i02 ■■■(f>nW (f)n+l • ■ ■ (f)N-l(f>N] , (4) 

where w = —a^ for r = +, w = Gz for r = 0, and w = cr+ for r = — . It 
is noted that when A= 1, $ represents the singlet state, and $n^\ 
and $n represent, respectively, the triplet states with 5'fotai ~ ^i 0' 
— 1. In the case of <C 1 the domain wall is expected to be trapped at 
the impurity bond, i.e., the iVth bond. We therefore use ■> 
and ^^j^^ as trial functions. For simplicity, we assume that 6i = 6 for all 
£'s. Then, these trial functions become equivalent to that used by Kennedy 
and Tasaki."*^ We have carried out a variational calculation with respect to 
9 for the energy expectation values. Up to the order of A^, the four trial 
functions lead to the same solution 6 = 6, where 6 is given by 



for D > Max(4, A) , 

iarccos(ff^ 

TT 
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6 = { l&TCCosi^^j for 4>L>>2A-4, (5) 

for Min(A, 2A - 4) > i:> . 



This solution in the region Z)>Max(4, A), that in the region 4>Z)>2A— 4, 
and that in the region Min(A, 2X — 4) > D correspond, respectively, to the 
large-D, Haldane, and Neel phases. (In the calculation for the Haldane 
phase we have neglected the factors cos^6 and cos^2^.) We denote $ and 
for ^ = 6* by l» and , respectively. It is noted that $ and for 



D = and A = 1 are equivalent to the valence-bond-solid^^ wave functions 
expressed by the Schwinger bosons. 

Let us concentrate our attention on the Haldane phase. For this phase, 
the variational results for Eg and Ai{K,) are given by 

-^-t^^ (6) 

Aoi^)c,{^^^^\n\^^^^)-E, 

f2{4-D){A-2X + D) ^ ~ _ 

= ^ { -k\^^Ao, (7) 

A,iK)c,{^P\n\^^^^)-E, 

= K 4^3^{2(8 ^2D- XD) - K{12 - 4A + D)} = kAi . (8) 

~ (t) 

When K = 0, as is expected, the four states $ and ^)^^ are degenerate 
in energy, which reflects the hidden Z2 x Z2 symmetry in the open chain. 
This degeneracy is resolved linearly with respect to in the case of /^>0, 
either < Ai(k) < Ao(k,) or < Ao(k) < Ai(k,) depending on whether 
D{4-\-2X — SK) is larger or smaller than 4(4— i^)(A — 1), and in the case of 
K<0, either > Ai > Aq (k) or > Aq > Ai (k) depending on whether 
D{4+2X-3K) is larger or smaller than 4(4-iC)(A-l). Thus, when A=l 
and D = 0, a transition of the ground state from a singlet to a triplet takes 
place at K, = with increasing k. This finding indicates that the impurity 
bond plays the role of the Z2 x Z2 symmetry breaking field and that the 
Z2 X Z2 symmetry is not a characteristic property of the bulk. 

The site dependence of the magnetic moments is calculated as 

($|5||$) = (I>SJV||I>5?)) = 0, (9) 
(l^^^l^lll*^^) = -{¥-^\SI\^^-^) = sin26^(cos^-i26^ + cos^-^26') . (10) 



The result given by eq. (10), which is essentially the same as that obtained 

by Kennedy and Tasaki,^) indicates that ($^'' 15*11 $^'') decays exponen- 
tially with the decay constant —(In | cos 29\)~^ as a function of the distance 
from the impurity bond. 

We briefly discuss the case in which the magnetic- field term Hmf = 

-^a;5'f„tai-^^'S'totai ('^'fotai = J2f=iS^) IS added to the Hamiltonian H. 
The energy difference AE (the energy measured from Eg) versus magnetic 
field Hx or Hz diagram can be constructed by solving the equation, 

[AoiK) - AE] [{Ai(/.) - A^;}' - -G^h'^{A,{k) - AE} = (11) 

with 

G = (^iv^^l^fotail^i?^) = (^iv"^l^fotail^i?^) = ^in 2^7(1 + cos22^~) . (12) 

We note here that ,,i|l>Sv^^) = (^l^fotail^J?^) = (^l^fotail^Sv ^) = 

($y'^|5fo^ai|^5v ^) = 0- '^^^ explicit results for the solution of eq. (11) will 
be discussed in a separate paper. ^) 

Next we deal with the case of k,<1, again concentrating our attention 
on the Haldane phase. In this case it is important to take into account the 
effect of the propagation of domain walls in the chain. When k,j^1, the 
scattering of propagating domain walls by the impurity bond produces a 
domain- wall bound state around it. To discuss this bound state, we employ 
Koster and Slater's method. (A more sophisticated treatment will be 
discussed elsewhere. ^^^) This method leads to the following result for the 
energy differences Ai{K,): 
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(13) 



where A^(l) and are constants, the latter being a quantity related to 
the matrix elements of /z-a^^i, and are determined phenomenologically below. 
The present Koster-Slater analysis also leads to the result that the energy 
difference between the bottom of the energy continuum for S'^otai ~ 
the ground state and that between the bottom of the energy continuum 
for 5'fQ^g^i = ±l and the ground state are given, respectively, by Ao(l) and 
Ai(l). Accordingly, these energy differences are independent of k,. 

We assume that eq. (13) also holds for ~ to determine Ai(l) and 
tti- Then, should be equal to Aj(l), since Aj(K) should vanish at k, = 0. 
Thus, we have 

A,(/s:) = A,(l){l-(l-/.)2} . (14) 

Furthermore, if we use the variational results given by eqs. (7) and (8), 
A^(l) is determined to be Ai(l) = Ai/2. 

On the basis of the results of the above Koster-Slater analysis, we 
discuss how the state at the bottom of the energy continuum changes with 
considering, for simplicity, the isotropic case of A = 1 and D = K = in 
which Aq(k,) = Ai(k,) = A(k,). At k = 1 the energy difference between the 
bottom of the continuum for the triplet state and that for the quintet state 
is considered to be equal to A(l). On the other hand, the binding energy 
of one trapped domain wall at «: = is given by A(l) [see eq. (14)], and the 
energy difference A(l) between the bottom of the continuum for the triplet 
state and the ground state is independent of k. Thus, when the value of k 
decreases, the state at the bottom of the continuum is expected to change 
at K = from the triplet to the quintet, since at k ~ the latter state is 
constructed by adding one trapped domain wall to the former state. 

Combining these considerations with the variational result that A(k) 
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is proportional to k, at k^O [see eqs. (7) and (8)] and also with the result 
given by eq. (14), we conjecture that the energy (measured from E^) versus 
K, diagram is that which is schematically given in Fig. 1. Here we have 
assumed that the energy difference between the bottom of the continuum 
and the ground state remains unchanged also for k<0. According to this 
diagram the Haldane gap is defined as the energy difference between the 
bottom of the energy continuum and the ground state for any k,. 

In order to numerically examine the analytical results as well as the 
conjecture discussed above, we have performed diagonalizations by Lanczos 
method for finite- chains in the isotropic case. In Fig. 2 we plot 
A(«:)/A(l) as a function of K'[={K-Kc)/{^-t^c)] for N = 12 and iV = 13, 
where k,c is the value of k, at which A{k,) vanishes for each N. The corre- 
sponding plot of eq. (14), where k,' = k,, is also given in Fig. 2. From this 
figure we see that the result for the k, dependence of A(«:) obtained by 
the above single domain-wall analysis is in rather good agreement with the 
numerical results. In Fig. 3 the energies (measured from Es) of low-lying 
excited states for 13 are plotted as a function of k,; the dotted line is 
for the lowest triplet state [i.e., the dotted line shows A(k)], the dashed 
line is for the second-lowest triplet state, and the dot-dashed line is for the 
lowest quintet state. This figure should be compared with Fig. 1. We plot 
in Fig. 4 the ^-dependence of the magnetic moments {SD of the S^^^^i = 1 
state for the lowest triplet and of the S^^^^^ = 2 state for the lowest quintet 
obtained for N = 13 and k, = (the open chain), where (• • •) stands for the 
expectation value with respect to a given state. This figure suggests a Neel- 
type configuration of the magnetic moments around the impurity bond for 
both states. It should be noted that the magnetic moments shown in Fig. 4 
are quite different from {Sf) = S^^^^g^^/N for k, = 1 (the periodic chain) and 
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that the values of the magnetic moments at the edges of the open chain 
are about 60 percent of the value 1 of the full moment. As may be seen 
from Fig. 4, the decay constant of the S^^^^i = 1 state for the lowest triplet 
is much longer than the variational result of —(In | cos 29 |)~^ ~0.9. In fact, 
Miyashita and Yamamoto"'"^-' have very recently shown that the N ^ oo 
value of the decay constant of this state is given by about 6. Details of the 
numerical calculation will be published in the near future. 

In conclusion, we emphasize that the Haldane gap should not be defined 
as the singlet-triplet excitation energy, but as the energy difference between 
the bottom of the energy continuum and the ground state. We propose, for 
the isotropic open chain, the energy (measured from E^) versus magnetic- 
field Hz diagram given by Fig. 5 (see Fig. 1 also). It is noted that the critical 
field H^, at which the finite magnetization per site appears, is equal to A(l) 
as in the case of the isotropic periodic chain. Furthermore, we expect that 
for samples doped by nonmagnetic impurities of a few or more percent, the 
magnetic response due to the Neel-type configuration of magnetic moment 
near the edges can be observed at a magnetic field even below H^. This 
may be an important subject for future experimental studies for examining 
our physical picture. 
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Figure Captions 

Fig. 1. Schematic energy (measured from Eg) versus k, diagram for the 
isotropic case obtained by analytical methods. The dotted line is for 
the lowest triplet state and the dashed (dot-dashed) line is for the 
bottom of the energy continuum consisting of triplet (quintet) states. 
The gray region shows the energy continuum. 

Fig. 2. PlotsofA(«:)/A(l) as a function of «:'[=(«:-A>:c)/(l-«:c)] foriV=12 
(open circles) and = 13 (open squares) in the isotropic case, where 
K,c is the value of k, at which A{k,) vanishes for each A^. The solid line 
is the corresponding plot of eq. (14), where k' = k,. 

Fig. 3. Energy (measured from Eg) versus k diagram for iV = 13 in the 
isotropic case obtained by numerical diagonalization. The dotted line 
is for the lowest triplet state, the dashed line is for the second-lowest 
triplet state, and the dot-dashed line is for the lowest quintet state. 

Fig. 4. Plots versus £ of the magnetic moments {S^) of the *S'^otai~-'^ state 
for the lowest triplet (full circles) and of the S'^otai ~ ^ state for the 
lowest quintet (open circles) obtained for AT" = 13 and ^^ = 0. 

Fig. 5. Energy (measured from E^) versus magnetic- field diagram pro- 
posed for the isotropic open chain. 
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